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ABSTRACT 



Aims. We verify whether the O/H abundances of galaxies can be derived from the equivalent width (EW) R23 instead of the extinction- 
corrected flux R23, and eventually propose a method of improving the reliability of this empirical method, which is often used for the 
non-flux calibrated spectra of galaxies. 

Methods. We select 37,173 star-forming galaxies from the Second Data Release of the Sloan Digital Sky Survey (SDSS-DR2), which 
offers a wide range of properties to test the EW method. 

Results. The EW-R23 method brings with it a significant bias: for the bulk of SDSS galaxies, it may affect the determination of 
log(0/H) by factors ranging from -0.2 to 0.1 dex and for some galaxies by factors ranging from -0.5 to 0.2 dex. We characterize this 
discrepancy (or bias) by a = (I[on]/lH^)/(EW[on]/EWH^), which is virtually independent of dust extinction, while tightly correlating 
with D n (4000), although at a lower significance, with (g - r) colors. 

Conclusions. The EW-R23 method cannot be used as a proxy for the extinction-corrected flux R23 method. From analytical third-order 
polynomial fits of a versus (g - r) colors, we have been able to correct the EW-^23 method. With this additional and easy correction, 
the EW-/?23 method provides O/H abundance values similar to those derived from the extinction-corrected flux R23 method with an 
accuracy of +0.1 dex for >92% of the SDSS galaxies. 

Key words. Galaxy: abundances - Galaxies: evolution - Galaxies: ISM - Galaxies: photometry - Galaxies: spiral - Galaxies: starburst 



1. Introduction 

Chemical abundance is a fundamental parameter for tracing the 
history of star formation and evolution of galaxies. Oxygen is 
the most commonly used metallicity indicator in the interstel- 
lar medium (ISM) by virtue of its high relative abundance and 
strong emission lines in the optical part of the spectrum (e.g., 
[O n]^3727 and [O m]/U4959, 5007). However, the "direct" 
method of estimating oxygen abundances from electron tem- 
perature (T e ) is generally only available for metal-poor galax- 
ies (12+log(0/H)<8.5), where the [O m]A4363 emission line 
can possibly be detected, which is needed for measuring T e by 
its ratio to [O ra]/U4959, 5007 (Pagel et al. 1992; Skillman & 
Kennicutt 1993). For metal-rich galaxies, the most commonly 
used are empirical strong-line methods, such as R23'. 



data sets from surveys. Then, equivalent widths (EWs) are be- 
ing used to replace the fluxes of their R23 values for estimating 
metallicities of the galaxies, i.e., from the EW^: 



EW R 23 



EW([0II]/13727) + EW([OIH]/U4959, 5007) 



(2) 



^23 = 



7([0II]/13727) + I([OIH]/U4959, 5007) 



(1) 



EW(H/3) 

This replacement was first checked by Kobulnicky & Phillips 
(2003; KP03 hereafter) on the basis of a sample of 243 
nearby galaxies. Consequently, this method has been used by 
some researchers to estimate metallicities of nearby and even 
intermediate-z galaxies (e.g. Kobulnicky et al. 2003; Kobulnicky 
& Kewley 2004; Lamareille et al. 2005a,b; Mouhcine et al. 
2006a,b etc.). 

However, it is known that there is a continuum (Fca) scale 
factor between the flux (F^) and EW (W^) values of the emission 
line: 



i.e. the flux ratios of [O 11] and [O m] to U/3 (Pagel et al. 1979; 
McGaugh 1991; Zaritsky et al. 1994; Kobulnicky et al. 1999; 
Tremonti et al. 2004 and the references therein). 

However, flux calibrations are frequently problematic in the 
current generation of wide-field galaxy surveys of multiobject 
spectrography, because of unfavorable observing conditions or 
instrumental effects such as a variation in system response over 
the field of view, nevertheless, one still expects to derive metal- 
licity properties of the star-forming galaxies detected in the large 



Fca 



(3) 
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This will naturally make us ask whether the continua underlying 
[O 11] and Hfi (and [O m]) are very similar or not. If not, Fca 
could be far from a unique constant, and then the metallicities 
derived from EW R23 could have large discrepancies from those 
derived from flux R23 . 

The SDSS provides a complete dataset and measurements up 
to several ten-thousand galaxies, making it a very good database 
for studying this question. In this paper, we selected 37,173 star- 
forming galaxies from the SDSS-DR2 to further check the reli- 
ability of using EW7?23 replacing the extinction-corrected flux 
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R23 to estimate the metallicities of galaxies. Moreover, this is a 
large homogeneous database observed by one single highly ef- 
ficient facility, which will minimize the effect of using various 
instruments. Some other characteristic parameters provided by 
the SDSS can also help for understanding the question further, 
e.g., Z) w (4000), g and r photometric magnitudes, etc. 

This paper is organized as follows. The sample selection cri- 
teria are described in Sect. 2. In Sect. 3, we check how big the 
difference is between the underlying continua of [O 11] and H/3, 
which is quantified as a parameter a (by 1 1 a). The emission-line 
quantities of the sample galaxies are analyzed in Sect.4, as well 
as the discrepancies between their EW7?23 and flux R23, and the 
discrepancies between the derived log(0/H) abundances from 
them. In Sect.5, we try to find the factors that mostly affected 
the a parameter, then to find the relations between them and a, 
hence to modify the EW R23 method, which includes the stellar 
population indicators D n (4000) and colors. In Sect.6, we discuss 
the boundary of the upper branch of 12+log(0/H) abundances 
from the \0gR23 calibration, then summarize and conclude this 
paper in Sect.7. 




a 8.5 9 9.5 



12+log(C5/H) (R,, 3 ) 

Fig. 1. Comparison between the oxygen abundances obtained by 
the MPA/JHU group (marked as 'SDSS') and those derived from 
R23 calibration of Tremonti et al. (2004) (marked as '^23')- The 
dashed lines mark the 0.1 dex discrepancies from the equal val- 
ues (solid line). 



2. Sample selection 

The data analyzed in this study were drawn from the SDSS-DR2 
(Abazajian et al. 2004). These galaxies are part of the SDSS 
"main" galaxy sample used for large-scale structure studies 
(Strauss et al. 2002). We selected a sample of star-forming galax- 
ies with metallicity estimates from the SDSS-DR2 database. The 
selection criteria are similar to those of Tremonti et al. (2004) 
and Liang et al. (2006). We summarize the criteria for selection 
as follows, and mark the selected number of the sample galaxies 
in parenthesis at the end of each step: 

(i) SDSS-DR2 (261,054), 14.5<r<17.77 mag (Petrosian magni- 
tude, 193,890 left); 

(ii) 12+log(0/H) S Dss > (50,385 left) (SDSS refers to the 
metallicity values provided by the MPA/JHU group, which were 
obtained on the basis of the 2001 Chariot & Longhetti model and 
Bayesian technique, see Tremonti et al. 2004); 

(iii) 0.04 < z < 0.25 (40,693 left) (this allows the spectral cov- 
erage from [O ii]/13727 to [S ii]/16731 for the sample galaxies); 

(iv) the emission lines [O 11], [O in], Ha, and [N 11] are de- 
tected and their fluxes and EWs are reasonable; fluxes of emis- 
sion lines Ha, and [N 11] 65 84 are detected at levels higher 
than 5cr (39,029 left); 

(v) 8.5<12+log(0/H)tf 23 <9.3 (38,932 left) (we adopt the formula 
of Tremonti et al. 2004 to convert R23 to the oxygen abundance 
12+log(0/H)# 23 , which is only suitable for metal-rich galaxies; 
therefore, we only select the galaxies with 12+log(0/H)>8.5 
here; 9.3 is almost the upper limit of the abundances in the sam- 
ples); 

(vi) the discrepancy between log(0/H)# 23 and log(0/H)sDSS is 
less than 0.1 dex (37,173 left finally), which removes some scat- 
tered data points from the R23 calibration, as shown by Fig.Q] 

Finally, we obtained a sample of 37,173 star-forming galax- 
ies from the SDSS-DR2 as our sample in this study. 

The criteria (i)-(iv) are almost the same as what was adopted 
in Tremonti et al. (2004), except the lower limit for redshift z 
was increased from 0.005 to 0.04 by following Kewley et al. 
(2005) to minimize the aperture effects of the SDSS. In Sect.O 
we discuss criterion (v) in particular, i.e. 12+log(0/H)~8.5 as 
the lower boundary of the upper branch of oxygen abundances 
from log/?23 calibration. 



3. Quantifying the difference between E\NR 2 3 and 
flux R23: determining the a parameter 

The MPA/JHU collaboration has put the measurements of 
emission-lines and some physical parameters for a large sample 
of SDSS galaxies on the MPA SDSS websiteQ (Kauffmann et al. 
2003; Brinchmann et al. 2004; Tremonti et al. 2004 etc.). These 
measurement values were obtained from the stellar-feature sub- 
tracted spectra with the spectral population synthesis code of 
Bruzual & Chariot (2003). 

Fluxes of the emission lines should be corrected for dust ex- 
tinction, so we estimate the attenuation of the sample galaxies 
using the B aimer line ratio Ha/Hfi and assuming case B recom- 
bination, with a density of 100 cm" 3 , a temperature of 10 4 K, 
and the intrinsic Ua/Uj3 ratio of 2.86 (Osterbrock 1989), follow- 
ing the relation of 

(T^U, = (^WlO^™. (4) 

1h(3 1HJ30 

Using the average interstellar extinction law given by 
Osterbrock(1989), we obtain f(Ua)- f(H0) = -0.37. For the 
data points with c < 0, c = is assumed since their intrinsic 
YLa/Uj3 may be lower than 2.86 if their electron temperature is 
high (Osterbrock 1989). 

Then, the extinction-corrected R23 parameter and its relation 
with the EW/?23 are: 

p Aoii] + Aonrj 

^23 - " 

_ Fiom 1Q c(f(roniw(Hi)) + [jom_ l0 c(Moiii])-f(m) 

Fh/3 Fuj3 

= W [0m F c ,[o m iQ 4 f(lom) _ f(m) 

, w [oin]^c,[oiii] 1Q c(monri)-f(H/3)) 

W\\pFc,wp 
a^[oii] + W[oiii] 

where a = (Fc 9 [om/F c ^)lO^^°^- f ^\ and 
(^Ctoin]/^,^)^^^ 0111 ^-^^ is about equal to 1 since 

1 http ://w w w. mpa- garching . mpg . de/S D S S/ 
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[O m] /L14959, 5007 and UJ3 are very close in wavelength. The 
expression /[oii],[oiii],h/3 refers to the dereddened, calibrated 
flux values of the corresponding lines; ^[oiiuoiiilh^ refers to 
the observed flux values; and W[oin,[om],H/8 re P resents me EW 
values of the related emission lines. The parameter c is the 
extinction coefficient. (Also see KP03) 

The a values of these SDSS galaxies can be calculated di- 
rectly from a = (F Ct [om/F c ^)l^ [0m) - f{H ^\ where F c ,[om 
and Fc,Hp can be estimated from their ratios of fluxes to EW val- 
ues, and c is the extinction coefficient. The derived a parameters 
of these sample galaxies show a median value of 0.85 and an av- 
erage value of 0.86 in a range from 0.1 to 2.6. KP03 estimate the 
typical value of the a parameters of their 243 sample galaxies to 
be a = 0.84 ± 0.3. However, KP03 still adopted a = 1 when they 
used the EW7?23 to estimate the oxygen abundances of galaxies. 

4. Analysis of emission-line quantities 

Several relations are analyzed in this section, including those of 
the emission-line quantities, e.g. EW(H/3), with the a parameter, 
with the discrepancies of EW R23 and flux R23 , and with the dis- 
crepancies of the log(0/H) abundances derived from these two 
R23 values. 

4.1. Relations between the line strengths and a parameters 

Figure [2^ shows the relations between the emission-line 
strengths EW(H£) and the line-ratio ^^-/^ (=l/c*) of the 
sample galaxies. The solid line marks the one-to-one correspon- 
dence. It seems that there is some correlation between the a pa- 
rameter and the emission-line strengths: from the galaxies with 
stronger emission lines, the ratios of log( ^°" ] / ^™ ) (=log(l/aO) 
change from -0.2 (mostly) to 0.5 monotonicly, though the scat- 
ters are large and show some exceptional points with log(l/a)<- 
0.2. This range is similar to what KP03 found for their sample 
galaxies. The line strengths of [O 11] and [O m], i.e., EW[oirj 
and EW[oni], also show similar trends to EW(H/3) in these 
kinds of relations, but with somewhat larger scatter. The ratio 
logCl^/l^ 1 ) ( =1 og(!M)) of all the sample galaxies show a 
median value of 0.069 dex (with a scatter of 0.060 dex) and a 
mean value of 0.077 dex (with a scatter of 0.104 dex). 



log(0/H)/? 23 as functions of the emission-line strengths EW(H/3). 
The general trend is that, from the galaxies with stronger line 
strengths, the differences in the two log(0/H) estimates change 
from -0.5 to 0.2 dex, with most of them ranging from -0.2 to 

0. 1 dex. KP03 do not present such a direct comparison for the 
O/H abundances. 

The large sample of SDSS star-forming galaxies that we use 
generally show that the EW7?23 will underestimate the oxygen 
abundances of the sample galaxies by a factor of 0.041 (the me- 
dian offset, with a scatter of 0.036 dex) or 0.054 dex (the mean 
offset, with a scatter of 0.078 dex). Here the log(0/H) abun- 
dances were obtained from the R23 calibration of Tremonti et 
al. (2004). We also adopted some other R23 calibration formulas, 

1. e. Kobulnicky et al. (1999, the analysis formulas for the mod- 
els of McGaugh 1991), Zaritsky et al. (1994), and Kobulnicky 
& Kewley (2004, the average of McGaugh 1991 and Kewley & 
Dopita 2002), to check these discrepancies, the results of which 
are quite similar. Moustakas & Kennicutt (2006) find similar dis- 
crepancies to ours by comparing the abundances derived from 
the extinction-corrected flux- and EW-/?23 of 12 nearby spiral 
galaxies. They find that the integrated abundances determined 
from the emission-lines systematically underestimated the char- 
acteristic abundances. The mean offset is -0.06+0.09 dex using 
the McGaugh (1991) calibration, or -0.1 1+0. 13 dex using the 
Pilyugin & Thuran (2005) calibration, and the corresponding 
median offsets are -0.04 and -0.09 dex, respectively. 

These results show that the oxygen abundances derived from 
the EW7?23 and the extinction-corrected flux-/?23 are not seri- 
ously different for these star-forming galaxies, generally less 
than 0.1 dex. However, the global discrepancy of them does 
show a wide distribution, from -0.5 to 0.2 dex, which means 
that this replacement of EW R23 to flux R23 for metallicity es- 
timates could cause different effects on the individual galaxies, 
and should be considered carefully. 

5. Modifying the EW /^-method 

It is interesting to further check the main factors that affect the 
discrepancy between the log(0/H) abundances derived from EW 
R23 and flux R23. We may then find some ways to modify the 
EW R23 method in order to obtain almost identical oxygen abun- 
dances to the flux R23 on the basis of this large set of SDSS 
galaxies. 



4.2. Discrepancies between EWR23 andR23 

When directly using the EW^3 to replace the flux R23, namely, 
a=l is adopted for EW%, the derived metallicities could 
have some discrepancies. Figure^ shows the discrepancies be- 
tween the EW7?23 and the extinction-corrected flux R23 as a 
function of line strengths EW(H/3). The general trend is that, 
from galaxies with stronger line strengths, the differences be- 
tween log(EW/?23) and log(7?23) increase from -0.15 (mostly) 
to 0.5 dex monotonicly. KP03 also find a similar trend for their 
sample galaxies. Our much larger sample shows this system- 
atic discrepancy more clearly. The median value of these dis- 
crepancies is about 0.061 dex (with a scatter of 0.050 dex), and 
the mean value of them is about 0.069 dex (with a scatter of 
0.086 dex). 

4.3. Discrepancies between (0/H) EWR23 and (0/H) Rl3 

One of the most important results of this work is given in 
Fig. [2fc, which shows the difference between log(0/H) £W Rl3 and 



5.1. Factors that affect the a parameter 

In Sect. 3, we point out that the a parameter is the ratio of the 
intrinsic continua underlying [O n]/13727 and Hfi and that it is 
related to the observed continua and the line extinction c (see 
Eq. ©). In this section, we discuss the dominant effect on the a 
parameter. 

Equation ([5]) shows a is a function of the dust extinction and 
stellar populations of the galaxies, so then we have 

a _ ^c,[oirj 1QC (f (roup- f (#/?)) 
Fcwp 
F° 

_ C,[OII] 10 (c-c*)(f([OIl])-f(Hj3))^ ^ 

where F^ [on] and F^ ^ are the dereddened continua underlying 
[O 11] and H/3, c the dust attenuation on emission-line, and c* 
characterizes the dust attenuation on the continuum. Here we 
assume c* and c follow the same reddening law. (Also see KP03) 
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Fig. 2. (a). Emission-line strengths EW(H/3) of the sample 
galaxies as a function of the ratios of equivalent width ratios 
to dereddened emission-line fluxes ratios for [O n]/H/3. (b). 
EW(H/3) as a function of the discrepancy between the quantity 
of EW#23 and the flux R 23 (extinction-corrected), (c). EW(H/3) 
as a function of the discrepancy between the two metallicity es- 
timates from EW R23 and flux R23 . 



As a rough estimate, we assume that the sample galaxies 
follow the same extinction law as the starburst galaxies studied 
by Calzetti et al. (1994), who found that the difference between 
the optical depths of the continua underlying the B aimer lines is 
about one-half of the difference between the optical depths of the 
Balmer emission lines (their Eq.26). Thus, we assume c* « 0.5c, 
and then obtain the following equation: 



_ c,[oii] 1 n o.5c(/([on])-/(#/?)) 



(7) 



where [ h]/^ch/3 characterizes the stellar populations of the 
galaxies, and c-term characterizes the dust extinction. 

It is easy to check the relation between a and the dust ex- 
tinction c since we have obtained both of them for the indi- 
vidual galaxies. However, we cannot obtain the intrinsic values 
°f [Oirj/^CH£ directly. Fortunately, the MPA/JHU group and 
the SDSS provide Z) w (4000) parameters and several photometric 
magnitudes for this large sample of galaxies, which can charac- 
terize the stellar populations of the galaxies. We check the rela- 
tions between a and dust extinction c, Z) w (4000), g — r colors for 
the sample galaxies in the next three sections. 



s 1 




Fig. 3. Relation of a parameter and dust extinction Ay for the 
sample galaxies. 



5.2. The relation between a and dust extinction 

The [O n]/l3727 is bluer in wavelength and is affected more 
strongly by dust extinction than thus the a parameter may 
be correlated with dust extinction Ay. Figure [3] shows the a pa- 
rameter as a function of dust extinction A v (=cRy/ 1 .47, Seaton 
(1979), R V =3A). The derived A v is from to 2.7, with the 
median and average values of 0.87 and 0.89, respectively. It 
shows that there is no clear correlation between a and Ay for 
these SDSS sample galaxies, and the linear least- square fit is 
a = 0.056Ay + 0.821 with a very slight slope. This means that 
the differences in the continua underlying [O 11] 3727 and Hfi are 
not affected much by dust extinction. The reasons may be that 
the related two lines are not very far away at wavelength and 
that the dust extinction coefficients of these SDSS star-forming 
galaxies are not so large. 



5.3. The corrected EWR23 method using D n (4000) index 

The break occurring at 4000A is the strongest discontinuity in 
the optical spectrum of a galaxy, and it arises because of the ac- 
cumulation of a large number of spectral lines in a narrow wave- 
length region. The main contribution to the opacity comes from 
ionized metals. In hot stars, the elements are multiply ionized 
and the opacity decreases, so the 4000- A break will be small for 
young stellar populations and large for old, metal-rich galaxies 
(Kauffmann et al. 2003). Kauffmann et al. (2003) adopt the nar- 
row definition from Balogh et al. (1999) and denote this index 
as £^(4000) for the SDSS galaxies. It is the ratio of the aver- 
age flux density F v in the bands 3850-3950 and 4000-4100 A, 
100 A narrower than the definition by Bruzual (1983). The pa- 
rameter Z) w (4000) is one of the main ones used by Kauffmann 
et al. (2003) to trace the stellar formation history of the SDSS 
galaxies, which shows a monotonic increase after the instanta- 
neous burst of star formation. 

We plot the relations of a against Z)„(4000) for the sample 
galaxies in Fig. [4^, which clearly shows a correlation. A third- 
order polynomial fit for this relation is obtained by fitting the 16 
median- value points in bins of 0.05 in £> w (4000) from 1 to 1.8, 
and given as: 



10.88 - 18.31 x + 11.18 x 2 - 2.34 x\ 



(8) 



with a standard error of 0.164, where x = D n (4000). This 
relation of a vs. £>„(4000) could be used to correct the 
EW R23 and then to obtain the consistent metallicities of galax- 
ies with the flux R23. We propose using ^23(EW) CO rrected= 
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Fig. 4. (a). Relation of the a and D w (4000) parameters for the 
sample galaxies. The big squares are the 16 median- value points 
in the bins of 0.05 in D„(4000) from 1.0 to 1.8. The solid line 
represents the third-order polynomial fit for these median-value 
points, and it is given as Eq. (8]). (b). Comparison between 
the metallicities derived from the ^-modified EW R23 by using 
D„(4000) following Eq. ® (marked as the subscript "modi- 
fied"). The diff refers to the difference between the oxygen 
abundances from the modified EW R23 and flux R23. (c). The 
direct comparison between the two oxygen abundances derived 
from the modified EW7?23 and flux R23 (extinction-corrected). 
In (b) and (c), the solid lines are the equal- value lines, and the 
dashed lines show the 0.1 dex discrepancy. 



axEW([0 n])/EW(H0) + EW([0 m])/EW(H0) to then estimate 
the metallicities of galaxies. 

Figure 0J) shows the consistency of the metallicity esti- 
mates from the corrected EW R23 and flux R23 after we apply 
the correction of D n (4000) for a. Figure [4k shows more di- 
rect comparisons for the two O/H estimates with such a cor- 
rection. Both of them show that the two derived abundances 
are very consistent. The median and mean discrepancies be- 
tween log(0/B) EW r 23 and log(0/H)# 23 now decrease to about 
-0.005 dex (with a scatter of 0.024 dex), and -0.010 dex (with 
a scatter of 0.054 dex), respectively. Then the corrected EW R23 
method provides log(0/H) abundances similar to those of the 
extinction-corrected flux R23 method within an accuracy of 
±0.1 dex for >94% of the galaxies. 

Unfortunately this correction is difficult to handle with non- 
calibrated spectra, simply because reliable measurements for 
Z)„(4000) require flux-calibrated spectroscopy. Since Z) w (4000) 
amplitude depends on the stellar population, age, and metallic- 



ity, it also correlates with colors. In the following we aim at gen- 
erating a correction usable by a large community, using a non- 
calibrated spectrum and one color as input (e.g. g - r for the 
SDSS data). 

5.4. The corrected EWR23 method using g-r color 

Colors can provide important information for the stellar popu- 
lations of the galaxies and can be obtained directly from the 
photometric observations. SDSS has made the u, g, r, i, z band 
magnitudes of the galaxies available publicly. Since the u mag- 
nitude has large uncertainty (20 per cent; Kauffmann et al. 2003) 
and the images in u- and z-band are relatively shallow, whereas 
/-band image may suffer from the 'red halo' effect (Michard 
2002; Wu et al. 2005), we use g-r color here to study such 
corrections for a. The g and r magnitudes can be converted to 
other band magnitudes, e.g. B, V, following some conversions, 
for example, Smith et al. (2002) and Jordi et al. (2006). 

Figure [5^ shows the correlation between the g - r colors (in 
Petrosian magnitudes) and a for the sample galaxies. The ba- 
sic trend shows that the redder galaxies have relatively lower a 
values, corresponding to the larger differences between the un- 
derlying continua of [O 11] and Hf3. The third-order polynomial 
fits for these correlations were obtained and given as 

a = 1.20 - 1.11 x+ 1.15 x 2 - 0.53 x\ (9) 

with a standard error of 0. 1 86, where x refers to g-r color. Then 
we propose using fl 23 (EW) corr ected= ax EW([0 n])/EW(H0) + 
EW([0 iii])/EW(H/3) to estimate the metallicities of galaxies. 

Figure [5j) shows the consistency of the metallicity estimates 
from the corrected EW R23 and flux R23 after we apply the cor- 
rection from g-r color for a. Figure [5k shows the compari- 
son between the two O/H estimates with this correction more 
directly. They show that the median and mean discrepancies be- 
tween log(0/H)Ew r 23 an d log(0/H)# 23 now decrease to about 
-0.004 dex (a scatter of 0.030 dex) and -0.012 dex (a scatter of 
0.062 dex), respectively. Then the two oxygen abundance esti- 
mates are almost identical to each other. Namely, the corrected 
EW R23 method provides log(0/H) abundances similar to those 
of the extinction-corrected flux R23 method within an accuracy 
of ±0.1 dex for >92% of the galaxies. 

5.5. Correction for the EWR23 method by a constant a=0.85 

The calculated median a value of this large sample of SDSS 
local star-forming galaxies is about 0.85 (the average value is 
-0.86). We may also use this constant ^-factor to modify the 
EW R23 . However, this constant correction is only useful for es- 
timating the global metallicity distribution of a large dataset, 
for example, from the database of a survey. As for the individ- 
ual galaxy, this may reversely enlarge the uncertainty for some 
of them, for example, the object with about EW(H/3)~10A(see 
Fig. [2k). Indeed, the correction on EW7?23 for the individual 
galaxies correlates tightly with their stellar populations and star 
formation histories. 

6. Discussions about the boundary of the upper 
branch of 12+log(0/H) 

In this work, we adopt 12+log(0/H)~8.5 as the low boundary of 
the upper branch of oxygen abundances from the log7?23 calibra- 
tion. The main reason is that the R23 calibration used (taken from 
Tremonti et al. 2004) is valid above this metallicity. In Sect. [6] 1, 
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Fig. 5. (a). Relation of the a and g — r colors for the sample 
galaxies. The big squares refer to the 10 median- value points in 
the bins of 0.1 in g - r from 0.1 to 1.1. The solid line represents 
the third-order polynomial fit for these median- value points, and 
it is given by Eq. ® . (b). Comparison between the metallicities 
derived from the a-modified EW R23 by using g - r colors fol- 
lowing Eq. © (marked as the subscript "modified"). The "cliff 
is as a function of g - r colors, (c). The direct comparison be- 
tween the two oxygen abundances from the Qf-modified EW R23 
and the flux R23 (extinction-corrected). The solid lines are the 
equal- value lines, and the dashed lines show the 0.1 dex discrep- 
ancy from the equal- value lines in (b) and (c). 



we present more observational data with O/H abundances de- 
rived from electron temperature T e , especially from the recent 
SDSS, and some photoionization models to further identify the 
reason we adopt 12+log(0/H)~8.5 as the boundary of the upper 
branch of oxygen abundances to compare the EW R23 and flux 
R23 methods. However, this boundary value is a bit higher than 
used by some other researchers, e.g., Pilyugin (2000, 2001a,b) 
and Pilyugin & Thuan (2005), who use 12+log(0/H)~8.2 on the 
basis of a sample of H 11 regions having their log(0/H) abun- 
dances estimated from T e . Therefore, we extend the boundary to 
compare the EW R23 and flux R23 methods and check whether 
our relations are valid or not in the region of 12+log(0/H)r e ~ 
8.2-8.5, which will be presented in Sect. [6)2. 

6.1. The observational data with (0/H) Te and the 
photoionization models 

In Fig. [6^ we plot the (0/H)r e -^23 relations for the sample galax- 
ies having T e -based O/H abundances, which are taken from 



Kniazev et al. (2004, 624 samples), Izotov et al. (2006, 409 sam- 
ples), and Yin et al. (2007, 695 samples). Figure 03 presents 
the predictions of photoionization models for the relations of 
12+log(0/H) vs. log#23 taken from Kewley & Dopita (2002) and 
Kobulnicky et al. (1999, K99), which was obtained by analyzing 
those ofMcGaugh(1991). 

Figure [6^ and Fig. [6j) show that both the large sample of 
observational data and photoionization models confirm that it is 
very difficult to derive a reliable O/H abundance from the R23 pa- 
rameter for the region of 12+log(0/H)<8.5 because of the large 
scatter of the data points there, the weak dependence of O/H on 
R23 (down to 12+log(0/H)<7.9), and the strong effects of ioniza- 
tion parameters. Yin et al. (2007) has used their Fig. 3b to show 
the large discrepancy, up to 0.4 dex, of the two sets of log(0/H) 
estimates derived from the upper branch and lower branch of 
R23 calibrations for the samples within 7.9<12+log(0/H)<8.4. 
Stasinska (2002) discusses the weak dependence of (O/H) on 
R23 in the transition region from nebular physics. 

Considering the discussions above, and also to avoid that 
some galaxies included may lie in the lower branch of the R23- 
(O/H) relation, we select the galaxies having 12+log(0/H)>8.5 
to check the EW R23 method in this study. However, it could 
be useful to check this and the validity of our relations in an 
extended upper branch range, 12+log(0/H)~8.2-8.5, as used in 
some studies. 

6.2. Checking the extended upper-branch range of 
12+log(0/H)~8.2-8.5 

We select a subsample from Fig.[6^ to check the situation in the 
extended range of upper branch, 12+log(0/H)~8.2-8.5. To be 
consistent with other parts of this work, we select this subsam- 
ple based on the SDSS-DR2 catalog. Finally, 37 independent ob- 
jects with metallicities of 8.2<12+log(0/H)r e <8.5 are selected 
by cross-correlating the DR2 catalog with the lists of Kniazev et 
al. (2004, from DR1), Izotov et al. (2006, from DR3), and Yin 
et al. (2007, from DR4). Then we use the EW and flux measure- 
ments of the related emission lines provided by the MPA/JHU 
group to estimate their EW R23 and flux R23, hence, the resulted 
abundances. 

These 37 galaxies show a discrepancy between log(EW7?23) 
and log(7?23) (as Fig. 2b) within a range of -0.2 to 0, which may 
result in overestimated log(0/H) abundances of about 0-0.2 dex 
by the EW7?23 as shown in Fig. 2c. Their £>„(4000) values are 
around 1.0, with g - r colors around 0.0 (ranging from -0.2 to 
0.3), which confirms that they are low-metallicity objects that 
will distribute in the left hand parts in Figs. 4a and 5a. Their a 
values are within 1-2, with the average value about 1.5. If we 
apply the a corrections for their EW R23, the correction factors 
will be about ar~ 1.4 and 1.2 by extrapolating Eqs. (8) and (9), 
respectively. 

If we do not consider the discrepancy among the differ- 
ent photoionization parameters and the weak dependence of 
O/H on R23 in this metallicity region (shown as Fig. 6b) and 
try to extrapolate the R23 calibration of Tremonti et al. (2004) 
down to these objects with 8.2<12+log(0/H)r e <8.5, then the 
open squares in Fig. 7a show that the R23 method will often 
overestimate the abundances of the galaxies and even under- 
estimate the abundances of some galaxies. The stars in Fig. 7 a 
show that the EW R23 method provides a higher log(0/H) abun- 
dance than the flux R23, generally about 0.2 dex, which is con- 
sistent with the discussions above and with Fig. 2c. Figure 7a 
also shows that the correction by a(g - r) will not improve 
the consistency for these objects much, and the correction by 
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log(R 23 ) 12+log(0/H) (flux R 23 ) 



Fig. 6. (a). The T e -based abundances and \0gR23 of the galaxies 
and H 11 regions taken from Yin et al. (2007, Yin07), Kniazev 
et al. (2004, Kniazev04), and Izotov et al. (2006, Izotov06). 
Associated with the calibration from SDSS-DR2 observations 
given by Tremonti et al. (2004; the thick solid line, marked 
by 'T04' ). (b). The predictions of photoionization models for 
the relations of 12+log(0/H) vs. log(7?23): the two dashed lines 
marked with y values are taken from Kobulnicky et al. (1999, 
K99), which are the analytical formulas for McGaugh (1991); 
the seven thin lines down to 12+log(0/H)~7.6 are taken from 
Kewley & Dopita (2002) with different ionization parameter q 
values. The thick solid line is the same as (a). Both (a) and (b) 
show that it is very difficult to derive the reliable O/H abundance 
from the R23 parameter for the cases of 12+log(0/H)<8.5. (See 
the online color versions of the plots for the discrepancies of the 
data and models.) 

a(D n (4000)) gives more consistent O/H abundances for the ob- 
jects with 12+log(0/H)/? 23 >8.2. We changed to the R23 calibra- 
tion of Zaritsky et al. (1994) (which is an average of three previ- 
ous calibrations) to obtain these comparisons again, and find the 
results are very similar to Fig. 7 a. 

We also used the R23 calibration of Pilyugin (2000) (their 
Eq.5) to do these comparisons, and present the results in Fig. 7b. 
It shows that mostly the flux R23 method provides consistent 
log(0/H) abundances with the T e method, except for some data 



Fig. 7. Comparisons between the log(0/H) abundances obtained 
from flux R23 and those from T e or from the corrected or uncor- 
rected EW R23 method: (a) using the R23 calibration of Tremonti 
et al. (2004, T04); (b) using the R23 calibration of Pilyugin 
(2000, POO). The different symbols represent those what ob- 
tained from the different methods marked by the labels in the 
top left corner. 



points in the left hand part of the data with underestimated 
log(0/H)# 23 and some reversed ones in the right hand section. 
It also shows that the corrected EW R23 method by <*(£>„ (4000)) 
provides more consistent abundances with the flux R23 for the 
objects with 12+log(0/H)^ 23 >8.1 than the uncorrected EW R23 
or corrected by a(g - r), which are similar to Fig. 7a. These 
differences among the EW R23 abundances presented in Fig. 7 
are not difficult to understand since the average a value of these 
low-metallicity objects is about 1.5, which is higher than the cor- 
rection relations provided, ~ 1.4 or 1.2. 

In short, in the lower metallicity range of 8.2< 
12+log(0/H)^ e <8.5, we would not recommend using our 
relations to correct the EW R23 method for the oxygen abun- 
dance calibrations, since there are several situations affecting 
the results, such as the big scatter of the data, the strong effects 
of ionization parameters, and the weak dependence of O/H on 
R23 in the range. 
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7. Summary and conclusion 

The goal of this paper is to check the reliability of using EW R23 
(= EW([QI ^^ ([Qm]) ) to replace the extinction-corrected flux R23 







_ /([OII])+/([OIII]) 

I(Hfi) 



) to estimate the metallicities of star-forming 



galaxies on the basis of a large sample (37,173) of galaxies 
with 12+log(0/H)tf 23 >8.5, selected from the SDSS-DR2. This 
replacement is often adopted when there are some problems 
dealing with proper flux calibrations for the spectral observa- 
tions. This large sample can provide some obvious statistical 
trends. 

The results show that the logarithm values of EW R23 and 
extinction-corrected flux-7?23 have a discrepancy from -0.4 to 
0.5 dex, with a median value of about +0.061 dex and a mean 
value about +0.069 dex. Thus, the discrepancies between the 
log(0/H) abundances obtained from EW R23 and those from the 
flux-7?23 range from -0.5 to 0.2 dex and have a median value of 
about -0.041 dex, a mean value of about -0.054 dex. These dis- 
crepancies are caused by the different continua (Fca) underlying 
the emission lines [O 11] and H/3 ([O m], as well). These differ- 
ences are characterized by the a parameter as (Fc,Hp)/(Fc,[oii]), 
which changes from 0.1 to 2.6, and by a median value of 0.85 
and a mean value of 0.86. 

Then we discuss the factor that affects this discrepancy 
mostly. Our large sample of data shows that the a parameter 
is almost independent of the dust extinction inside the galaxies 
and depends closely on stellar populations of the galaxies, which 
can be quantified by the Z)„(4000) parameters and colors of the 
galaxies. Third-order polynomial fits have been obtained for the 
observed relations of a versus Z)„(4000) and a versus g-r colors 
for the sample galaxies, which can be used to modify the EW R23 
method. After applying this correction by D n (4000), the median 
and mean discrepancies between log(0/H)EWR 23 an d log(0/H)# 23 
decrease to about -0.005 dex and -0.010 dex, respectively. After 
applying this correction by g-r colors, the median and mean dis- 
crepancies between log(0/H)^^ 23 and log(0/H)# 23 decrease to 
about -0.004 dex and -0.012 dex, respectively. The two derived 
sets of log(0/H) abundances are almost identical now. 

In summary, when there are problems with flux calibrations 
of the spectra, the EW R23 could be used roughly to replace the 
extinction-corrected flux R23 to estimate the metallicities of star- 
forming galaxies. The discrepancy caused by this replacement 
can be from -0.2 to 0.1 dex generally. This factor is consistent 
with those found by KP03 and Moustakas & Kennicutt (2006). 
However, this discrepancy could be large for the different indi- 
vidual galaxies, from -0.5 to 0.2, if the underlying continua of 
[O 11] and UJ3 ([O m]) are quite different. The D„(4000) parame- 
ters and colors of the galaxies are very useful for correcting the 
EW R23 method, which will then greatly decrease the discrepan- 
cies and result in an almost identical oxygen abundance to the 
flux-7?23. 

Nevertheless, the modified EW7?23 still suffers from the 
drawback of the "double- value" of R23 for the O/H estimates, 
and some other strong-line ratios are also useful for esti- 
mating the metallicities of galaxies then, such as [N n]/Ha, 
([O m]/H/?)/([N n]/HoO, [N n]/[S 11] and [O in]/H/3, etc. (Kewley 
& Dopita 2002; Pettini & Pagel 2004; Liang et al. 2006; Yin et 
al. 2007). 
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